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Abstract 



We suggest a new mass generation mechanism for gauge fields. The 
quantum field theory constructed in this paper is nonabelian, gauge in- 
variant and asymptotically free. 
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It is well known that the naive mass generation mechanism can not be applied 
in case of the nonabelian Yang-Mills theory. We recall that if F^ v — d^A^—d^A^ 
is the strength tensor of the abelian U(l) gauge field A^, fi = 0, 1, 2, 3, and S is 
the usual Maxwell action, 



then one can easily construct the action Sf^ generating the corresponding mas- 
sive theory (see, for instance, [8], Sect. 3.2.3), 



Here and thereafter we use the standard convention about summations and low- 
ering tensor indexes with the help of the standard metric g^ v of the Minkowski 
space, g o = 1, 9u = -1 for i = 1, 2, 3, and = for i ^ j. 




(1) 




+ ^-A^\d i x. 



(2) 
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In contrast to (1) action (2) is not invariant under gauge transformations 

A tl ^A„ + d^ip, (3) 

where ip is any scalar function on the Minkowski space. However, if to 7^ the 
equations of motion generated by action (2), 

0"F„„ + m 2 A u = 0, (4) 

imply the Lorentz gauge fixing condition 

= 0. (5) 

Condition (5) is obtained by taking the divergence of equation (4). 

Since action (2) is quadratic in A^ the field theory associated to (2) can be 
easily quantized, the key point in the quantization procedure being the gauge 
fixing condition (5) which automatically appears in the massive case. 

If G is a compact simple Lie group the construction of the massive quantum 
vector field theory outlined above can not be carried over in case of the Yang- 
Mills field associated to G. Formally the action similar to (2) can be written 
down. But the new action containing selfinteraction terms turns out to be non- 
renormalizablc (see [8], Sect. 12.5.2). The main obstruction to renormalizability 
in the nonabelian case is nonlinearity of the corresponding gauge action. This 
leads, in turn, to the fact that the gauge fixing condition similar to (5) does not 
automatically appear in the nonabelian case. 

Usually one overcomes this difficulty by using spontaneous symmetry break- 
ing and the Higgs mechanism [5, 6]. But as a byproduct of that mechanism one 
gets extra unwanted boson fields. 

In order to avoid the Higgs mechanism for constructing a massive nonabelian 
gauge quantum theory we start by slightly modifying action (2) in such a way 
that the new action S m for the massive field is invariant under gauge transfor- 
mations (3). We introduce the new action S m by the following formula 

S m = J (-J V" + ^(PAMPA)^ d 4 x. (6) 

Here 

(PA)„ = Ap — d li n~ 1 ff'A v (7) 

is the orthogonal projection of the field A^ onto the kernel of the divergence 
operator d ■ A = d^A^, with respect to the scalar product < •, • > of vector 
fields, 

<A,B>= j A t ,(x)B f "(x)d 4 x, 

and □ = d^d^ is the d'Alambert operator. In formula (7) it is assumed that 
an inverse to the d'Alambert operator acting on the image of □ is fixed. This 
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can be achieved, for instance, by choosing boundary conditions for large xq = 
t. The corresponding boundary conditions should be also imposed on gauge 
transformations tp in formula (3) to insure that each <p belongs to the image of 
□ and the identity 

D- 1 Dip = ip (8) 

holds. The particular choice of the operator D _1 is actually not important. In 
this paper we shall assume that all inverse d' Alambert operators are of Feynman 
type. 

Observe that (PA)^ is the gauge transform (3) of for = —\3r 1 d v A v . 
By definition the quantity (PA)^ is invariant under gauge transformations 
(3). Indeed, if A^ A^ + d^<p then 

(PA)„ ^A^ + - d ll n- 1 ff'{A v + hip) = (PA)„. 

(PA)n is actually the gauge invariant part of A^. Moreover, (PA)^ satisfies the 
Lorentz condition (5), 

d^(PA)„ = d»A„ - d >1 d fJ n- 1 d v A v = 0. (9) 

Note that (PA)^ linearly depends on A^ since the gauge action is linear in the 
abelian case. 

Since (P A)^ is invariant under gauge transformations action (6) is also gauge 
invariant. 

Observe that the lagrangian density associated to (6) is not local. From 
the first sight this is a serious obstruction for using action (6) in quantum field 
theory. However, if 

WAp = 

one has by definition (7) 

(pa),, =a„- d^n-'d-Av - v 

Therefore after imposing gauge condition (5) action (6) coincides with action 
(2) and, hence, generates the same equations of motion as (2) modulo gauge 
transformations. We conclude that action (6) can be quantized similarly to (2), 
and the corresponding Green functions for local gauge invariant observables 
constructed with the help of (6) satisfy the locality condition. 

From the above discussion we infer that actions (2) and (6) are physically 
equivalent and generate the same quantum field theory In this sense action (6) 
is a trivial generalization of (2). But a remarkable property of the former one is 
that it has a gauge invariant renormalizable counterpart in the nonabelian case! 

Formula (6) is not yet the most suitable for generalizing to the nonabelian 
case since it contains the operator P which is related to the explicit form of 
the gauge action (3). In order to exclude the operator P from formula (6) we 
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observe that since (PA)^ is a gauge transform of A^, and the strength tensor 
is gauge invariant, one also has 

F„ v = d»A v - d v A^ = d^(PA) u - d„{PA)^. (10) 

Now using (9), (10) and integrating by parts we obtain the following formula 

-J ^{U-^F^F^^x = J (PA)^PAYd i x. 

The last relation allows to rewrite action (6) in an alternative form, 

S m = ~\j ( V + m^n^F^F^) d*x. (11) 

The linearity of the gauge action in the abelian case is actually not crucial for 
formula (11) since it does not explicitly contain the operator P. 

Now it is easy to write down a nonabelian counterpart of action (11). First 
we fix the notation as in [8]. Let G be a compact simple Lie group, g its Lie 
algebra with the commutator denoted by [•,•]. We fix a nondegenerate invariant 
under the adjoint action scalar product on g denoted by tr (for instance, one 
can take the trace of the composition of the elements of g acting in the adjoint 
representation). Let i a , a = 1, . . . , dimg be a linear basis of g normalized in 
such a way that tr{t a t b ) = -\5 ah . 

We denote by A^ the g-valued gauge field (connection on the Minkowski 
space), 

A — A a t a 

Sip — Si^L . 

Let be the associated covariant derivative, 

where g is a coupling constant, and F^ v the strength tensor (curvature) of A^, 

F^ = d„A„ - d v A tl - g[A„, A v \. 

We shall also need a covariant d'Alambert operator Da associated to the 
gauge field A^, 

n A = d^d*. 

The covariant d'Alambert operator can be applied to any tensor field defined 
on the Minkowski space and taking values in a representation space of the Lie 
algebra g, the g- valued gauge field A^ acts on the tensor field according to that 
representation. Note that the operator Oa is scalar, i.e. it does not change 
types of tensors. 

Finally recall that the gauge group of G- valued functions g(x) defined on 
the Minkowski space acts on the gauge field A^ by 

An^-id^g-i+gApg- 1 . (12) 
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The corresponding transformation laws for the covariant derivative and the 
strength tensor are 

D^gD^g- 1 , (13) 
F^^gF^g- 1 . (14) 

Formula (13) implies that the covariant d'Alambert operator is transformed 
under gauge action (12) as follows 

DA^gDAg- 1 . (15) 

In the last formula we assume that the gauge group acts on tensor fields ac- 
cording to the representation of the group G induced by that of the Lie algebra 
0- 

From formulas (14) and (15) it follows that the natural nonabelian analogue 
of action (11) invariant under gauge action (12) is 

Sg = ljtv (F^ + m 2 (a^F^)F^) d 4 x, (16) 

where in the expression (D~^ 1 F tll/ ) we assume that F M „ is in the adjoint repre- 
sentation of q and the scalar operator D^ 1 acts on F M „ componentwise. The 
first term in (16) is simply the action of the Yang-Mills field. 

Now a new problem immediately appears. Due to nonlinearity of the gauge 
action in the nonabelian case it is impossible to make the lagrangian density 
associated to action (16) local by imposing a gauge fixing condition. However, 
the appearance of nonlocal quantities in quantum field theory is not a new fact. 
One should recall, first of all, the Faddeev-Popov determinant [3] that can be 
made local using additional anticommuting ghost fields. In this paper we are 
going to apply a similar trick to action (16). 

Similarly to the case of pure Yang-Mills field (see [8], Sect. 12.2; [4], Ch. 3, 
§3) let us first formally write down the expression for the generating functional 
Z(J) of the Green functions associated to action (16) via a Feynman path 
integral, 

Z(J) = J V(A)V( V )V(fj) exp{ 4 J H\F, V F^ + ^-(D^F^F^ (17) 

+2rj(Mf]) - 2J^A ll ) + T(A)}d 4 x}, 

where J M is the source for taking values in the Lie algebra g, J 7 (A) is a 
gauge fixing term, M. the corresponding Faddeev-Popov operator, 7j, r\ are the 
anticommuting scalar ghost fields taking values in the adjoint representation of 
(Faddeev-Popov ghosts). Here and later we always assume that the measure 
V in Feynman path integrals is properly normalized. 

The key observation is that the nonlocal part of action (16) is the quadratic 
form of the operator D^ 1 . Therefore to make the expression in the exponent in 
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formula (17) explicitly local it is natural to use a formula for Gaussian integrals 
relating the exponent of the quadratic form of the operator and the exponent 
of the quadratic form of the inverse operator. For operator K acting on scalar 
fields that formula takes the following form ([4], Ch. 2, §6) 



ex 



ex - 



J (x-VMM = (is) 

= det{K)* J V(iP)exp{i J (-^(Kip)ip + ipip)d 4 x}. 

Since in formula (17) the operator D^ 1 acts on the skew-symmetric (2,0)- 
type tensor field F^, in order to apply a formula similar to (18) to the expression 

p{ , J tr(^(D^ F^F^d'x) (19) 

one has to introduce a new skew-symmetric (2,0)-type tensor field (bosonic 
ghost) with values in the adjoint representation of the Lie algebra g. 

Using the ghost field $^ and formula (18) we can rewrite expression (19) 
in an alternative form, 

exp{i J tr(^(D^ F^)F^)d 4 x} = (20) 

= (detn A \ QHg) )i J P($)exp{* J tr( - i(D A $^)$^ + ^$^ F ^)d 4 x}, 

where the operator Da is assumed to act componentwise in the space f2 2 (g) of 
skew-symmetric (2,0)-type tensors with values in the adjoint representation of 0. 
The numerical factors in the r.h.s. of the last formula appear since the Feynman 
path integral in (20) is only taken with respect to the linearly independent 
components $ M „, [i < v of the skew-symmetric tensor 

The r.h.s. of (20) still contains a nonlocal term, (detCU|n2( g ))3 . But since 
SI 2 (g) is the direct sum of six copies of the space fi°(g) of scalar fields with values 
in the adjoint representation of g and the operator CU acts componentwise in 
fP(g) the determinant detCU|n2( ) is equal to (detCU|fiO( )) 6 , and 

(detCU|fi2 (fl) )5 = (deta4ho (B) ) 3 . 

Therefore recalling the Faddeev-Popov trick (see [4], Ch. 3, §3; [8], Sect. 12.2.2) 
and introducing six scalar anticommuting ghost fields r\i,r\i, i = 1,2,3 taking 
values in the adjoint representation of q we can express the square root of the 
determinant, (detCU|Q2( B )) 2 ; in the following form 

(detCU|n*( fl ))* = (detDA|oo( fl )) 3 = (21) 
3 3 
= J n^)^)cxp{z J 2tr(J2v l (D A r ll ))d i x}. 



»=i " 1=1 
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Now substituting (20) and (21) into into (17) we derive an expression for the 
generating function Z( J) that only contains local interaction terms and suitable 
for developing perturbation theory, 

3 

Z(J) = / V(A)V(*)V(r,)V(Tf) nPfaiPfri) x 
J i=i 

x exp |i y [tr Qi^F"" - i(D A $^)^ + (22) 

+2^(M,7) + 2^^(n A% ) - 2J"A„^ + ^(^)]d 4 x| . 

Formula (22) looks a little bit complicated. But it is equivalent to the com- 
pact expression (17). Note that when deriving (22) starting from (17) we only 
used tricks with Gaussian Feynman path integrals. From (17) it also follows that 
when m — the generating function Z( J) is reduced to that of the quantized 
Yang-Mills theory. 

Now we can discuss renormalization of the lagrangian in formula (22). In 
order to do that we have to first fix a gauge, i.e. we have to choose a gauge 
fixing term T(A) and specify the Faddeev- Popov operator M. 

We shall see in a moment that for our purposes the most convenient choice 
of T(A) is 

F{A) = tr(0M„0M„ + m 2 9"i,(D- 1 9 , '4)). (23) 

Since the r.h.s. of (23) only depends on the divergence d^A^ of A^ the corre- 
sponding Faddeev-Popov operator has the same form as in case of the Lorentz 
gauge (see [4], Ch. 3, §3), 

M - d^ZV (24) 

Note that gauge fixing condition (23) contains a nonlocal term. But since 
the nonlocal term only depends on the longitudinal component of the gauge 
field A^ the nonlocal part in (23) does not generate extra nonlocal terms in 
the expressions for gauge invariant physical observables (see [4], Ch. 3, §3 for 
discussion of a similar situation in case of pure Yang- Mills field). 

Finally combining (22), (23) and (24) we arrive at the following formula for 
Z(J) 

3 

Z{J) - / V(A)V(^)V(r])V(fj) HviriVfc) x 

x exp ji J tr Qi^i^ - J(0 A Vr + 7? V " 2J ^+ (25) 
+277(^^7?) + 2j2m(OAVi) + 0M„0M„ + m 2 & t A li (n- 1 ff'A v ) \ d 4 x\ . 
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By simple dimensional counting the quantum field theory with generating 
function (25) is renormalizable. Moreover, using dimensional regularization and 
the BRST technique one can show, similarly to the case of the Yang-Mills field 
(see [8] ) , that gauge invariance is preserved by the renormalization. The details 
of those proofs and of other calculations that we are going to discuss can be 
found in [14]. 

The particular choice (23) of the gauge fixing term is convenient since in this 
case the free propagator D ab v of the gauge field has a very simple form, 

n ab _ A cab 9fJ,u 

^ D + m 2 ' 

Therefore the gauge field becomes massive. However, the theory may also 
contain some massless particles corresponding to the ghost fields. 

The coupling constant renormalization calculated at the one-loop order us- 
ing dimensional regularization coincides with that for the pure Yang-Mills field 

(see [8]), 

11 „2 S 



go^l-^C-j, (26) 

where g is the bare coupling constant, e — > when the regularization is re- 
moved, fi is an arbitrary mass scale, and the constant C is defined by 



CS ab = ^2 Ca d G'l 



cd 
b 



c,d 



with the help of the structure constants C° h of the Lie algebra 55, [t a , t b ] — C^ b t c . 

It is easy to see that the coupling constant renormalization in the theory 
with generating function (25) calculated within the dimentional regularization 
framework coincides with the coupling constant renormalization for the pure 
Yang-Mills theory to all orders of perturbation theory. Indeed, by the results of 
[7] in the case of dimensional regularization the coupling constant renormaliza- 
tion is independent of the mass for dimensional reasons. Therefore the coupling 
constant renormalization in the theory with generating function (25) is the same 
as in the massless case, i.e. in case of the Yang-Mills field. 

One can also calculate the mass renormalization at the one-loop order, 

mo = m^(l-4-2^C-), (27) 



16tt 2 6 e, 

where mo is the bare mass. 

The wave function renormalization at the one-loop order has the same form 
as in case of the Yang-Mills field, 

z *- 1 + iS?i4 (28) 
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From (26), (27) and (28) one can get the corresponding renormalization 
group coefficients at the one-loop order (see [13], Sect. 4.5), 

«e>=*i = -i&ir c ' < 29 > 

<91nm g 2 11 „ 

^)=^ = -^ 2l rC, (30) 

1 8lnZ A g 2 5 

7d(g)= 2^^^ = -16^3 C - (31) 

The Fourier transform T(pi, . . . ,p n ,m, g, /u) of any one-particle irreducible 
Green function with engineering dimension c? is transformed under rescaling of 
the momenta as follows (see [13], Sect. 4.6) 

T(sp 1 ,...,sp n ,m,g,iJ,) = (32) 

-n J — 7d(g(s')) 

where s > is a scale for the momenta. In equation (32) g(s) and m(s) are the 
solutions to the following equations 

3^=/?(g(s)), g(l) = g; (33) 



9s 
dm(s) 



m(s)(7 rn (g(s)) - 1), m(l) = m. (34) 



Applying standard renormalization group arguments (see, for instance, [8], 
[13]) we have g(s) — ► when s — ► oo since the beta function /3(g) is negative in 
(33), and the theory with generating function (25) is asymptotically free. 

Similarly, m(s) — ► when s — > oo since in that case g(s) — > as we have 
just observed, and hence by (30) 7 m (g(s)) — 1 is negative in (34) for large s. 
Therefore the mass can be neglected for large momenta. 

In conclusion we remark that the mass term in action (16) can be regarded 
as the result of a renormalization of the Yang-Mills action. Such possibility was 
discussed in [2]. However, the mass term can not be observed perturbatively 
at large momenta. We also note that there should be no serious problems 
with generalizing to our theory the program of rigorous construction of Green 
functions for the quantized Yang-Mills field outlined in [12]. 

One can also add to action (16) a fermionic part coupled to the gauge field 

in the minimal way. 

Acknowledgements. When this work was completed the author has 
received information that the gauge invariant mass term that appears in formula 
(11) in the abelian case was considered in [10], and action (16) in the nonabelian 
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case was recently introduced and studied in [1]. Paper [1] also contains an 
extensive list of references on dynamical mass generation in the Yang-Mills 
theory. 

Another action describing massive gauge fields was proposed in [11]. 

Prof. R. Jackiw also communicated to me that in case of three dimensions 
the mass generation mechanism similar to the one introduced in this paper was 
suggested in [9]. The author would like to thank him, David Dudal, Amitabha 
Lahiri and Jean-Luc Jacquot for pointing out the references quoted above. 
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